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Abstract 

The generalized Henon-Heiles Hamiltonian H = 1/2 (Pj- + Py + ciX"^ + 
C2Y^) + aXY"^ — 6X^/3 with an additional nonpolynomial term fiY^'^ is 
known to be Liouville integrable for three sets of values of (b/a, ci, C2). It has 
been previously integrated by genus two theta functions only in one of these 
cases. Defining the separating variables of the Hamilton-Jacobi equations, 
we succeed here, in the two other cases, to integrate the equations of motion 
with hyperelliptic functions. 



I. Introduction 



The generalization of the Henon-Heiles Hamiltonian defined by 

X' = Px, Y' = Py, (2) 
X" = -aY^ + bX^ -ciX, (3) 

Y" = -2aXY - C2Y + -^At, ^ arbitrary, (4) 

16a^ Y-^ 



is integrable in the sense of Liouville in three casesJM namely 

b 

- = -1, Ci = C2, (5) 

a 

- = —6, ci,C2 arbitrary, (6) 

- = -16, ci = 16c2, (7) 
a 

and is equivalent to the reduction ,^ = a; — ct of three fifth order soliton equa- 
tions,§ respectively the Sawada-Kotera (SK), KdVs and Kaup-Kupershmidt 
(KK) equations, by applying the translation 

SK: u = X + ^, c = ciC2, (8) 

KdVs: u = X + ^^^^, c= -^(-3c?-16ciC2 + 48cl), (9) 
2Ua IIJ 

KK: u = X + —, c = ciC2. (10) 
2a 

The reduced partial differential equations (PDE's) are respectively 

SK: + lOann" + — a V - cu + 4aS - ^ = 0, (11) 
3 3a 

KdVs: ti^^) + 20aun" + 40a^u^ + lOan'^ - cu 

+iaE + T-^(ci + 4c2)(c? - 12ciC2 + 16ci) = 0, (12) 
lUUa 

KK: ti(^) + 40ann" + — a\^ + 30au'^ -cu + AaE - ^ = 0, (13) 

3 3a 

in which E is the constant value of the Hamiltonian H. Therefore we will 
further use these names for refering to the respective integrable cases of the 
Hamiltonian. 
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The problem which we address is to find the separating variables and to 
explicitly integrate. 

The general solution for the KdVs case has been obtained in terms of 
hyperelliptic functionJl^ by separation of the variables of the Hamilton- 
Jacobi equation in parabolic coordinates. 

For fi = 0. the equations of motion for the SK and KK cases have been 
integratedErllSI in terms of elliptic functions. 

In this paper we give the general solution of the equations of motion for 
the SK and KK cases in terms of hyperelliptic functions. This is achieved 
by separation of the variables of the Hamilton-Jacobi equation, where the 
canonical transformation between SK and KK0 plays an important role. 

In section we consider the SK case with /i = as a starting point to 
treat the general case {fi arbitrary). In section ID, we recover the canonical 
transformation between SK and KK, starting from the factorization of the 
scattering operator associated with these equations and using the link with 
the Fordy-Gibbons equation.llHI In section [V, we give the general solution for 
KK, using the canonical transformation previously obtained. In section 0, 
we use the fact that the canonical transformation is invertible for obtaining 
the general solution of the SKcase. In both cases, in the limit /x ^ 0, we 
recover the previous results J3't3 In section VI, we compare our method of 
integration, starting from the Hamiltonian system, with the method used 
by Cosgroveta for integrating the fourth order ODE's (|Tl]) and 



II. Particular case: Sawada-Kotera for jj, = 

This is the simplest case and it will be our starting point to define the 
separation of variables of the Hamilton-Jacobi equation for KK. 

Consider the Hamiltonian in the SK case derived from ([l|) by setting 
a = ^ and defining ,U = X + C2,V = Y,c = ciC2 

H ^ Ki,o = liP^ + P^) + ^UV^ + - |C/, (14) 
which is a constant of motion of the equations 



U" = -UV^ + U^) + 



c 



2^' ' ^ ^ ' 2' (15) 
V" = -UV, 

where U,V and the derivatives U' = Pu,V = Py are functions of the 
independent variable ^ (' = ^). 
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This system possesses a second constant of motionll 



i^2,o = -2Pc/Py _ [/V - ^ + cV, (16) 

which is in involution with the Hamiltonian, i.e. {i^i,0) -^2,0} = 0. 
The separation of variables is defined as 



Qi = [/ + Q^ = U-V, 

Pi = i(Pc/ + Pv), P2 = - Pv), 

and in those new variables the expressions of Ki^ and ^2,0 are 



(17) 



i^i,o = P? + P| + ^iQl + Ql) - + Q2), (18) 

i^2,o = 2(P| - ) + ^(Q^ - Ql) - |(Q2 - Qi), (19) 



such that the equations of motion become 



I Q'l' = -Wi + cQi-K2,o + 2Ki^o, 
\ Q2=-IQI + cQ2 + K2,0 + 2Ki,o- 

They possess the general solution 

Ql = -Up U - 6, tI' -134(2^^1,0 - K2,o)) = -12pi(e), 

Q2 = -i2p (e - ^2, ii, -133(2^1,0 + ^2,0)) = -i2p2(e), 



(20) 



(21) 



in terms of the Weierstrass elliptic function p{x — xq, 02,93), solution of the 
first order differential equation 

p'2 = 4p3 _ _ (22) 
such that the solutions become 

/ [/ = -6(pi(0 + P2(6), ..o^ 

\ y = -6(pi(0-P2(0)- ^ ^ 

III. Canonical transformation between SK and KK. 

Let us consider, for /i arbitrary, the equations and the constants of mo- 
tion associated with the SK and KK Hamiltonianslli'0 

SK: U' = Pu, V' = Pv, (24) 
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U" = -\{V' + U')-'-, V" = -UV+^, (25) 

Ki = 1{P^ + P'v) + \UV' + -1^ + ^, (26) 

Kl = Kl, + ^fiU + fi^, (27) 

KK: ~ ^' ci = 16c2, c = ciC2, u = X + 2c2, w = y, (28) 

u' = pu, v'=pv, (29) 

It" = -^v^ - Au^ + c, v" = --uv + 4-' (30) 

h = ^{pl+pl) + ^uv^ + ^u''-cu + -^, (31) 

11 1 c 

^2,0 =pt- j^v^ - ^2^^^"^ + ^^^^^ ~ gPuPt^t^^ + Y^w^, (32) 

fci = 4o + f« + 2/^§ + ^- (33) 

The reason why the expressions A;2,o and k2 are defined by their square will 
appear soon. 

The two nonlinear partial differential equations SK and KK 

SK: Ut + {a,^^^ + 5UU^^ + ^U% = 0, (34) 

20 o 2 

KK: ut + {uxxxx + Wuuxx + y « + 30ti^)s = 0, (35) 



whose reductions ^ = x — ct are (O) and ([1^), respectively obtained from 
the systems (p5|) ~(p^) and (^0|)-(|3l|) by elimination of the variables V and 
V, possess a Lax pair with a third order scattering problem (L^lJ = Xip). 
The scattering operators can be factorized in the following way 

SK: L = d^, + Udx = idx-w)idx+w)dx, (36) 
KK: L = d^, + 2udx + ux = idx + w){d,)idx-w), (37) 

such that the solutions of the PDE's are related through a Backlund trans- 
formation 

SK: U = Wx- (38) 
KK: u = -Wx - -w'^. (39) 
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with the solution w of the Fordy Gibbons equationli 



XX - 5w'^Wxx - 5wwl + W^)x = 0. 



(40) 



The reduction ^ = j; — ct of this equation can be solved for w(^), either 
by eliminating U and w' between (|3^) and the equations of motion (|2^) 



1 



w 



2 y2 



(41) 



or by eliminating u and w' between (39) and the equations of motion (pQ) 



w 



-a V 



such that, defining = —fj, and 

r = 6{VK2,o + XPu), 

+2,X^uv^ + ISX'^v'^pI + UX^vp^ + SA"^), 
the canonical transformation is given by0ill 



Pv ^ X 



V 2^2 



Pu 



y3 



(3P{^ + 3[/y2p^ _ p^y3 



3A 
2^6 



A2 



1 



U 



Pv 



4y2 

-6 



2Pv + -]^-X 



V 



V 



Pv X 

9 

V 



n 

4y8 ■ 



1 



(12p^ + Guv'^py - v^pu) 



3A 

+-^{2uv^ + 12^2^2 + i2Awp„ + 4A2), 



1 



•y + -)VU+X 



(42) 

(43) 
(44) 

(45) 



(46) 
(47) 



(48) 
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IV. General solution of the Kaup-Kupershmidt case 



Starting from the separation of variables (|T^ and the canonical transfor- 
mation (^7|)-(^) in the case = 0, we consider the transformation defined 
byS on the basis of Painleve analysis 



Qi 

pi 

P2 



-6- 

-6 
1 



«2,0 



pI + hfi 



u, 



2v'^ 
1 

2^' 



(49) 



{I2pl + 6uv^pv - v^pu + I2pvk2,t 



This inverts to 



u 



P2 -Pi 



.q2-qi 



1 



{qi + 92), 



qi - q2 



Pu = 24 



Pi -P2 



qi - q2 



+ 2{pi -P2, 



qi + q2 , ^Piq2 -P2qi 



(50) 



+ 2 



pI 



12fe2,0 



{P2 -piY 



qi - q2 



qi - q2 



(qi - q2f 

Taking account that k2fi is no more a constant of motion, this change of 
variables will appear be useful to find the general solution for KK. In those 
new variables, the Hamiltonian system (|3l| ) becomes 



H 



^2,0 



9l 



kl=pl+pl + ^{qf + ql) - {{qi + ^2) + (51) 



'^{pl-pi) + \{ql-ql)-l{q2-qi), 



2pi + 



fj- {qi - q2)pi 



6 



"'2,0 



/ „ fj- (gi - q2)P2 
^2 = - 6 ^1 



^2,0 



Therefore, defining f{qi,Pi) = pf + ^qf - f( 
Jacobi equation is separated 



(52) 
(53) 

(54) 

1,2), the Hamilton- 



h[f{qi,Pi) - f{q2,P2)) = f{qi,Pi) - f{q2,P2) + j^iqi - 92), (55) 
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K = f^ (i = l,2)- 
We can express the second constant of motion k"^ in two equivalent ways 



or k2 = --q2 + {k2,0 - rrr—, , 

3 12 k2fi 



(56) 
(57) 



so that the ehmination of /x(gi — (/2)/^2,o between (|5lD, (^) and (|5lD, ( |57| ) 
yields 



^2 



fl3 

2 yi 



^qi + i-4pi-^ + cqi+2k,f, 



or A:2 = -f'?2 + (4p^ + |-cg2-2fci)2. 



(58) 
(59) 



The elimination of between (pq) and (|5^), and of p2 between (|59|) and 
(II) yields 



91 - 92 



9i 



9i 



= ^/2fci - f + cq2 + ^jkl + I 1 - 



91 - 92 



'(60) 



•(61) 



In the case /i = 0, the differential equations for qi and q2 are separated and 
their solution is expressed in terms of the Weierstrass elliptic function 

91,0 = -12p {i-ii, ^, -Yi^(2A:i,o - A:2,o)) ^ -12pi(0, (62) 
(?2,o = -12p - 6, ^, -Yi^(2^i,o + ^2,o)) = -12p2(0, (63) 
so that the solution for (|0D in the case = is0 

3 /p;(e)-p'2(0^' 



2 Vpi(e)-P2(0 

^2,0 

p2(0 - Pl(0 



+ 6(pi(o + P2(e)), 



(64) 
(65) 



In the case / 0, let us introduce the new variables 



S2 



'3^ +92, 



which transform the equations (|60|) and (61) into 



4 



i(.f-3M)3 + e(.2_3M) 



1 



^2 > 



/2fci-i(.2_3M)3 + c(si-3 



Si - S2 
1 



Defining 



2A;i - - I - 3^ 



+ c 



3 VS1-S2 



3^' 



(66) 

(67) 
(68) 

(69) 



the system (p7|)-(|6q) can be solved by inversion of the hyperelliptic integrals 



ds 



00 ^/P{s) 
1 sds 



00 



ds 



00 /P(i) 
2 sds 



00 



(70) 
(71) 



which define si and S2 as multivalued functions of 

The general solution of the equations of motion (|30|) in the case 7^ is 

u = -l{sl + sl) + hl-l('-l±^)\ (72) 
2 /i 2 Vsi + S2/ 

.2 = (73) 

Sl + S2 

As they are rational symmetric combinations of si and S2 and their deriva- 
tives, u and v'^ are single-valued functions of ^. 

In the variables qi , 52 this solution is expressed as 



2k^ 
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+ cqi - JB + ^gi 



2fci 



9i 



+ Cg2 + 



V 



gi - q2 



{qi + Q2), 



^2 + I9l + V^2 + 1Q2 

qi-q2 



J 

(74) 
(75) 
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which clearly goes to (|6^-(|6^) in the limit — > 0. 

V. General solution of the Sawada-Kotera case 

We start from the general solution (^), (^) for KK and apply the canon- 
ical transformation ( |47| ) to obtain the general solution for the SK Hamilto- 
nian system 



U = V^(s'i + 4) + sl + siS2 + sl- -Kl (76) 
= -2V^{si + S2){sis[ + S24) + 2(si + S2f I si + sl- ^ I • (77) 



2/x 

Let us check that the limit of this solution when /i — > is (p3|). The 
expression (^) for U can also be written as 

2 /i 2 V-si + S2. 

+ + • ^^^^ 



_3 / s[ + 4 
Since in the limit ^ — > 



one has 



Si + S2 = 0{^), 

S'l+S'2 ^ pm-p'2io 

S1+S2 Pl(0-P2(0 



limU=l{Q, + Q2). (79) 



Next, for the expansions 

2(S1 + S2)' ( + ^2 - ^ j = 4(^1 - Q2f + Oifi), /i ^ 0, 

{si + S2)(sis'i + S2S2) = /" ^ 0, 

provide the limit of V"^ in ([77|). Therefore, 

limV^ = \iQi-Q2f. (80) 
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VI. Comparison with the results of Cosgrove 

C. Cosgrovelll recently integrated the ODE's ( [Tl| ) and (|13|) with hyper- 
eUiptic functions, using the postmultipher method. 

To compare the two different ways of integration, let us recall the vo- 
cabulary introduced in Painleve analysis of nonlinear differential equations, 
making the distinction between fixed and movable constants. A constant 
is called fixed if it appears explicitly in the differential equation, while it is 
movable if it is a constant of integration and therefore depends on the initial 
data. 

In the Hamiltonian formalism, described by the system ([2^)-(p3D and 
(P^ -(P^), ki (resp. Ki) and /c2 (resp. K2) are movable constants, while is 
fixed (it appears in the equations of motion) . 

In Cosgrove's paper, the first integrals of the fourth order equation he 
obtained in formulas (4.3)-(4.4) and (5.6)-(5.7), which are therefore mov- 
able constants, correspond in the Hamiltonian formalism to k2 (resp. K2) 
and respectively introduced as movable and fixed constants. In order to 
integrate the resulting second order sixth degree differential equation and 
transform it in a coupled system of first order equations for applying the 
postmultipher method, Cosgrove defined "suitable" auxiliary variables cho- 
sen in a "subjective" way. Here the link between the canonical variables and 
its expressions (4.5)-(4.6) and (5.3)-(5.4) can be clearly established. Those 
expressions, which are "hidden" variables for the fourth order differential 
equation, are nothing else than the canonical variables v"^ (resp. V'^) and Py 
(resp. Py) explicitly defined in the Hamiltonian formalism. 

VII. Conclusion 

We have proven that the three integrable cases (SK, KdVs, KK) of the 
Henon-Heiles Hamiltonian can be integrated in terms of hyperelliptic func- 
tions. 

We will take advantage of the method developed here to integrate other 
Hamiltonian systems with two degrees of freedom and additional nonpoly- 
nomial terms. 
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